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ABSTRACT 
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fT") ■ Context. Compact binary supersoft X-ray sources (CBSS) are explained as being associated with hydrostatic nuclear burning on 

£N| ' the surface of a white dwarf with high accretion rate. This high mass transfer rate has been suggested to be caused by dynamical 

instability, expected when the donor star is more massive than the accreting object. When the orbital period is smaller than ~6 
hours, this mechanism does not work and the CBSS with such periods are believed to be fed by a distinct mechanism: the wind- 
i-C 1 driven accretion. Such a mechanism has been proposed to explain the properties of objects like SMC 13, T Pyx and V617 Sgr. One 

observational property that offers a critical test for discriminating between the above two possibilities is the orbital period change. 
Aims. As systems with wind-driven accretion evolve with increasing periods, some of them may reach quite long orbital periods. The 
above critical test may, therefore, also be applied to orbital periods longer than 6 hours. CAL 87 is an eclipsing system in the LMC 
with an orbital period of 10.6 hours that could provide the opportunity for testing the hypothesis of the system being powered by 
wind-driven accretion. 

Methods. We obtained eclipse timings for this system and show that its orbital period increases with a rate of P/P = +7.2(±1.3) x 10 s 
years. 

Results. Contrary to the common belief, we conclude that CAL 87 is the first confirmed case of a wind-driven CBSS with an orbital 
period longer than 6 hours. The system is probably an evolved object that had an initial secondary mass of M2, = 0.63M© but is 
currently reduced to about M 2 = 0.34Mq. 

Conclusions. We discuss evidence that other CBSS, like CAL 83 and V Sge stars, like WX Cen, are probably also wind-driven 
systems. This may in fact be the rule, and systems with inverted mass ratio, the exception. 

CO ' 

• ■ Key words, binaries: close - Stars: winds, outflows - binaries: eclipsing - Stars: individual: CAL 87 - supernovae: general - 

C ' Techniques: photometric 

O 

, 1 . Introduction hours and the mass donor star is certainly less massive than the 

■ white dwarf. In this case we have a mass ratio that is normal for 

'O ■ Com P act blnar y supersoft X-ray sources (CBSS) are a group of catacl ic variableSj and the dynamical instability does not oc- 
X ; objects, first discovered in the Magellanic Clouds, with unusual cur Hqw C£m we Mn ^ ired w fa accretion rate? Van 
g , properties. Their large (-Eddington) luminosity is mostly radi- Xeeseli & Ki (m ® have proposed that high mass trans- 

a ted in the supersoft X-ray spectral range of 20 to 80 eV. Van den fe rate couW bg driyen from a j irradiated low-mass 

Heuyel et al GUI have shown convincingly that this emission donQr staf Such an lanation was pr0 p 0S ed for SMC 1 3 (van 
is originated from hydrostatic nuclear burning on the surface of Xeeseli et al rjggg and for T p (Kni et al 2000), and 
a white dwarf, In order for this to happen a high accretion rate qquU ^ bg Ued tQ GQ Mug ^ ^ j t . nova hase 

(~ 1O- 7 M yr- 1 ) is required, such as may be realized in systems (Diaz & m9 ^ & Steiner [i5^. A critical test that 

with inverted mass ratios, in the sense that the mass donor is the couW bg Ued tQ ^ tWQ of sms h tQ measure ^ QI _ 

more massive of the binary components Mass transfer is then ^ riod ch a . ^ (M h mass) d mical instability mass 
dynamically unstable and occurs on the Kelvin-Helmholtz time- (DIMT) model edicts afl iod decrease while 

scale, with the orbital period decreasing with time For mass ±e (1qw mass) wind . driven mass transfer (W DMT) model pre- 
donor stars with masses in the range 1 .3 to 2.5 M Q , this configu- dicts a . d increase 
ration will provide the required accretion rate of ~ 10 MQyr . 

This framework has been adopt ed in t he literature since it was This critical t est wa s applied to the V Sagittae star V617 
proposed ( van de n Heuvel et al. [19921 see also Kahabka & van Sgr (Steiner et al. 12006). V Sge stars are cons idered the galac- 
den Heuvel [T997l for a review and references). tic counterpart of the CBSS (Steiner & Diaz 119981 1. All three 

The short period CBSS 1E0035.4-7230, also known as SMC wel1 known stars of this class have eclipses, and this provides 
13 (Schmidtke et al. [1996b . does not fit into this classical pic- a g° od opportunity for measuring variations of the orbital pe- 
ture of supersoft X-ray sources as being driven by dynamically riod - V Sge has an orbital period of 12.3 hours and a decreasing 
unstable mass transfer. This binary has an orbital period of 4. 1 orbital P eriod (Patterson et al. |1998|); this is consistent with the 

DIMT scenario. V617 Sgr has an orbital period of 4.98 hours 

* Based on observations made at the SOAR Telescope. (Steiner et al. 1 19991 1 and, therefore, is below the critical period 
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of 6 hours for the DIMT model (Deutschmann 1 1 9981 King et al. 
1200 11 1. In fact, an increasing orbital period has recently been re- 
ported (Steiner et al. 120061 1. confirming the WDMT scenario. In 
such a case, the orbital period increases with time, as the sec- 
ondary "evaporates", and may reach values as high as 30 hours 
(van Teeseling & King [l998l . 

One may, therefore, ask whether the other known CBSS stars 
with orbital periods longer than 6 hours are in fact DIMT cases 
or advanced stages of WDMT? Of course, each case would have 
to be tested in order to verify this dichotomy. Perhaps the most 
promising case to be studied in order to answer this question is 
CAL 87, a well known binary in the LMC (Cowley et al. [19901 
Hutchings et al. 119981 1. Accurate timings have been measured 
by Cowley et al. (|1990|) and, 7.1 years later, by Alcock et al. 
( 119971 1, allowing a precise determination of the orbital period. As 
9.8 years have elapsed since the last timing, with an additional 
determination one could make a second accurate measurement 
of the period, allowing assessment of any change. In the present 
paper we report timings for the eclipses of CAL 87 and show 
that the star has an orbital period that increases with time. 

2. Observations and data analysis 

CAL 87 was observed with the SOAR Optical Imager (SOI) on 
six nights in October and November 2005, during the Brazilian 
time of the early science phase of the SOAR Telescope. The data 
consist of time series of images obtained with the mosaic of two 
EEV 2048 x 4096 thinned and back-illuminated CCDs, through 
the Johnson V filter, and with individual exposure times of 60 
seconds. Each series spanned about 1 hour and was centered on 
the expected time of the eclipses. The data reduction - bias and 
flatfield corrections - was performed with IRAF0 standard rou- 
tines. Differential photometry of the images was executed with 
the aid of the DAOPHOT II package. 

From the six observed light curves, two were discarded be- 
cause of the presence of flares that disturbed the measurements 
of the eclipse timings. The remaining four eclipse timings were 
measured individually, as well as combined in a single set folded 
in phase with the ephemeris from Alcock et al. J 1997b . The av- 
erage timing obtained from our data is presented in Tabl^TJ in 
addition to the eclipse timings found in the literature. 




P(hr) 



Fig. 1. P/P as a function of the orbital period. Top - WDMT 
systems: the solid curve represents the expected loci for systems 
at the beginning of the wind-driven evolution. The dashed lines 
(from left to right) are the evolutionary tracks for systems with 
initial secondary masses between 0.1 and 0.6 Mq, in steps of 
0.1 M Q . The square shows the WDMT binary CAL 87 - see 
Table [5] Bottom - DIMT systems: the solid curves (from left to 
right) represent the expected loci for systems with primary star 
masses of 0.55, 1.0 and 1.4 Mq. The triangle shows the DIMT 
binary V Sge. 



P 2 = 0.442 677 66(5) d with respect to cycle 5884; this corre- 
sponds to the mean cycle E 2 =9915. The difference between P 2 
and Pi is 9 to 10 times the uncertainty in the determination of 
those periods. From the difference between P 2 and Pi and the 
elapsed time P.(E 2 - Ei) we obtain P = + 1.7(+0.3) x 10~ 10 , the 
uncertainty of which was estimated with the usual error propa- 
gation method from the errors of Pi and P 2 . As expected, the 
new timing detemination allowed the measurement of the or- 
bital period change: the period increases with a time-scale of 
P/P = +7.2(+1.3) x 10 6 years. 



Table 1. Eclipse timings of CAL 87. 



Obs. minimum 


E 


Reference 


(HJD) 


(cycles) 




2 447 506.8021(2) 





Cowley etal.|1990| 


2 450 111.5144(3) 


5884 


Alcock etal.|l997| 


2 453 680.8244(3) 


13947 


this work 



The ephemeris of Alcock et al. (11997b provides a period of 
Pi = 0.442 677 14(6) d\ this period was found from observa- 
tions spanning a total interval of 5884 cycles, with the mean cy- 
cle of Ei = 2942. From our timings we can derive a period of 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



3. Assisted stellar suicide - a terminal case? 

The first conclusion is that, contrary to the canonical view that 
considers CBSS as having massive (1-2 Mq) secondary stars, 
CAL 87 clearly has a mass ratio which is typical for cataclysmic 
variables, i.e., M 2 < M\. This has been anticipated by Cowley 
et al. (119901 119981 > who found, given the small amplitude of the 
radial velocity of emission lines, that M 2 - 0.4Mq and M\ > 
AMq. For this reason those authors explored the idea that such 
a massive primary could be a black hole. We now know that the 
mass accreting star is probably a white dwarf: its mass, derived 
from fitting the luminosity-temperature diagram with theoretical 
calculations, is M\ = 1.35M© (Starrfield et al. l2004l >. The ob- 
served supersoft X-ray luminosity, Lx = 4x 10 36 erg/s (Starrfield 
et al. 120041 1. implies an accretion rate of M acc = 10~ 8 Mq yr . 
This is additional information which can be used to solve the 
following equations (van Teeseling & King 119981 see also eqs. 
3, 4 and 5 in Steiner et al. l2006l >: 

jq_^M^P 
1 + g M, P 
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and 



1 +q) 



1.2 x 10- 6 M, 



where 

3y0 2 - q - 2 
: (l+g)(2-3 9 )" 



(2) 



(3) 



Here /?2 is the ratio of the specific angular momentum of the 
wind to that of the secondary star. We have also assumed that 
M2 < 0.6 (see van Teeseling & King 1998 for a detailed discus- 
sion). 

The right-hand side of equations ([TJ and (|2j are measured 
quantities, while from the left-hand side two independent param- 
eters can be determined: q = 0.25 and fii - 0.89. This implies 
that the mass of the secondary is M2 = 0.34Mq. 

One should explore the range of parameters that could 
be accommodated with the observations. The observed super- 
soft X-ray luminosity means that M acc cannot be smaller than 
10~ 8 M© yr~ x ; therefore the minimum values for ySo and q are 
0.89 and 0.25, respectively. A larger value for M acc could be 
possible if part of the X-rays are absorbed, for instance, due to 
the high inclination. However, equation (01 dictates that q < 2/3. 
For this reason M acc < 1.9 x 10~ 8 M© yr~ l , but M 2 could be as 
high as M2 = 0.9 M©. For all the range of possible parameters, 
we always obtain P2 = 0.89. 

Is this value for /?2 satisfactory? We adopt the idea that the 
wind from the secondary star must be blown from its heated face. 
We will make the approximation that the wind's average specific 
angular momentum is that of the point which is at R2/3 from 
the center of the secondary star, and ~ 2/3 of R2 from the inner 
Lagrangian point. This is approximately 



#, ~ 1 - (0.13 + 0.07 log q){\ + q) 



(4) 



Here we have assumed that the radius of the secondary is 
given by the well known Paczynski's relation Rn/a = 0.38 + 
0.2 log q (Paczynski 119711 ), where a is the separation between 
the two stars. For a mass ratio of q = 0.25, we obtain fio = 0.89, 
in agreement with that determined above. For higher values of 
the mass ratio, fio could be as low as P2 = 0.80 (for q = 2/3). 
The numbers derived above are, therefore, consistent; they also 
imply that, for a binary with an orbital period of 10.6 hours, the 
secondary is quite evolved. From van Teeseling & King ( 119981 1 
we have 



Pi * 3 ^P orb M2 hr 

and 
Mr, - Pi/9 



(5) 



(6) 



where P, and M2, are the orbital period (in hours) and the sec- 
ondary mass (in solar masses) at the beginning of the wind- 
driven evolution, respectively. 

This means that the initial orbital period was f,- = 5.7 hr 
and the initial mass of the secondary was about Mu ~ 0.63M©. 
In other words, 0.29M© have already been evaporated or trans- 
ferred to the primary star. Considering M\ = 1.35M© and the 
accretion rate of about 10~ 8 M© yr~ l , CAL 87 may explode as 
a type la supernova in 10 7 yr. The parameters of CAL 87 are 
presented in Table [2] 

Figure [TJpresents the PjP versus P relation for WDMT and 
DIMT systems. The curve for WDMT systems at the beginning 



of wind-driven evolution was calculated from Eqs. ([TJ, (f3j and 
(|6j and from the values of P2 and q derived above for CAL 87. 
The parameter PjP is highly sensitive to the value of /?2, so the 
curve plotted is not extensive to other WDMT systems. Each 
WDMT evolutionary track was derived from Eq. ([TJ and Eq. (|5j. 
To construct the curves of the DIMT systems, we assumed a 
mass-radius relation for main sequence stars for the secondary, 
and the expression for the radius of the Roche lobe of the sec- 
ondary from Paczynski (|19711 l, and derived the orbital period 
through the determination of the separation between the binary 
components. The values of PjP were obtained assuming mass 
and angular momentum conservation in binary systems 



PjP 



M\M 2 



3M 2 (M 2 - Mi) 



(7) 



and adopting the thermal timescale mass transfer for CBSS 
(Kahabka & van den Heuvel [T997j : 



M 2 = 



3 x 10 7 M 2 ' 



(8) 



Figure[TJillustrates the distinct behavior of the period change 
in the two regimes, WDMT and DIMT. 

For a mass of M2 = 0.34M© for the secondary, the pre- 
dicted radial velocity semi-amplitude is K\ = 68 km s _1 . This 
is consistent with the amplitude of the radial velocity as deter- 
mined by the absorption lines, but not consistent with the emis- 
sion lines amplitude (Hutchings et al. 1 19981 Cowley et al. 1 1 9981 >. 
Absorption lines also have the right phasing but large uncertain- 
ties. The emission lines have smaller amplitudes; Hell 4686 A 
has^i = 30± 8 km s " 1 (Hutchings etal.[jl98j and 40+ 12km s -1 
(Cowley et al. [T990j, w hile O VI presents K r = 35 ± 26 km s _l 
(Hutchings et al. |19981 >, but all are out of phase by about A<p = 
0.14. As all the lines have different amplitudes and phasing, 
we conclude that they result from a combination of absorp- 
tion lines superposed on emission, distorting their profiles; non- 
orbital motions in the disc, and wind emission, may also signifi- 
cantly affect the determination of the radial velocity amplitudes. 
Further studies are necessary to answer this question. 



Table 2. Physical parameters and orbital elements of CAL 87. 



Parameter 


Value 


Reference 


Mi 


1.35M© 


Starrneld et al.|2004| 


M 7 


0.34M© 


this work 


H 


0.25 


this work 


i 


78° 


Mever-Hofmeister et al.|1997| 


Porb 


10.6 hr 


Alcock etal. [19971 


M 2i 


0.63M© 


this work 


Pi 


5.7 hr 


this work 


Pi 


0.89 


this work 


PjP 


+7.2(±1.3)x 10 6 years 


this work 


T 2 


< 1200 K 


this work 



As the system presents deep eclipses, one could argue that, 
at mid-eclipse, the spectrum of the secondary star could reveal 
its evolutionary status. However, this may not be the case as its 
temperature is likely to be very low. The expected temperature is 
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Table 3. CBSS with measured P/P. 



Star 


P orb 


P/P 




Reference 




(hr) 


(xlO f 


years) 




V617 Sgr 


4.97 


+ 1.1 


+ 0.2 


Steiner etal.|2006| 


CAL 87 


10.62 


+7.2 


+ 1.3 


this work 


VSge 


12.34 


-5.4 


+ 0.2 


Patterson etal.119981 










Szaszetal.|2006| 



T2 < 1200 K. This value is derived assuming that the luminosity 
is smaller than the value for the main sequence of this mass, and 
considering the radius derived from the current system parame- 
ters. 

The inclination of the system was determined by Meyer- 
Hofmeister et al. (1997) assuming a secondary more massive 
than the primary. They found i = 78°. In the current scenario the 
inclination is higher, being i > 75° in order for the white dwarf 
to be eclipsed, but more likely 2 79° in order to produce the 
deep eclipse. The inclination cannot be much larger than this, 
otherwise we would not detect significant X-ray emission out- 
side the eclipse because of strong attenuation by the accretion 
disc. 

Table 4. CBSS with estimated mass ratios. 



Star 


Type 


Port (hr) 


q = M 2 /M l 


Refs. a 


RX J0537.7-7034 


CBSS 


3.48 


0.5 


1 


1E0035.4-7230 


CBSS 


4.10 


0.3 


2 


V617 Sgr 


VSge 


4.97 


0.4 


3,4 


WXCen 


VSge 


9.98 


0.35 


5,6 


CAL 87 


CBSS 


10.62 


0.25 


7,8 


VSge 


VSge 


12.34 


3.8 


9,10 


QR And 


CBSS 


15.84 


0.4 


11, 12 


CAL 83 


CBSS 


24.98 


0.4 


13 



References: 1- Greiner et al. 120001 2- van Teeseling et al. [l998l 3- 
Steiner et al. [T999l 4- Ste iner et al. l2006l 5- Diaz& Steiner [l995l 6- 
Oliveira & S teiner 120041 7- Hutchi ngs et al. [19981 8- this w ork; 9- 
Herbig et al. [T965l 10- Smak et al. ljOoTI 11-M cGrath et al. 1207)11 
12- Deufel et al. [T999l 13- Cowley et al. [T998l 



Are there other evolved WDMT systems among the long pe- 
riod CBSS, or is CAL 87 an exception? Cowley et al. ( TIWSl ) 
suggested that CAL 83 has a low-mass secondary, given its low 
amplitude radial velocity. Similarly, WX Cen shows indications 
that it has a low mass ratio (Diaz & Steiner |1995t Oliveira & 
Steiner l2004l l. From the known systems that belong to the class 
of CBSS and V Sge stars, six are short period (< 6 hr) systems 
(GQ Mus, T Pyx, RX J0439, RX J0537, IE 0035, and V617 Sgr) 
and are therefore WDMT systems. Among the long period ob- 
jects, four have indications of being WDMT systems (WX Cen, 
CAL 87, QR And and CAL 83, see Table|4]), and one has no in- 
dication to date of its mass (RX J0513), while MR Vel probably 
has a subgiant secondary (Schmidtke et al. |2000t . 

The only clear case for a DIMT to date is V Sge (Herbig 
et al. 119651 Patterson et al. 1998). This supports the argument 



that wind-driven mass transfer (with ten objects) may be the rule 
and dynamical instability mass transfer (with one object), the 
exception. Another implication is that the CBSS are an old pop- 
ulation, instead of a population of intermediate age as previously 
believed. 
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